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. metry breaking required to obtain the correct Z mass is quantified. To 
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fields is presented. This model has additional triplets in the low energy 
theory which get a mass of 0(500) GeV from a coupling to a singlet. 
The improvement in fine tuning is quantified and the spectrum in this 
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complete model is presented and some phenomenological constraints 
are discussed. 
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1 Introduction 



One of the outstanding problems of particle physics is the origin of elec- 
troweak symmetry breaking (EWSB). In the Standard Model (SM), this is 
achieved by one Higgs doublet which acquires a vacuum expectation value 
(vev) due to a negative mass squared which is put in by hand. The SM has 
the well known gauge hierarchy problem [Q. It is known that supersymme- 
try (SUSY) stabilises the hierachy between the weak scale and some other 
high scale without any fine tuning if the masses of the superpartners are less 
than few TeV P, The Minimal Supersymmetric Standard Model (MSSM) 
is considered as a low energy effective theory in which the soft SUSY breaking 
terms at the correct scale are put in by hand. This raises the question : what 
is the origin of these soft mass terms, i.e., how is SUSY broken ? If SUSY 
is broken spontaneously at tree level in the MSSM, then there is a colored 
scalar lighter than the up or down quarks [^]. So, the superpartners have to 
acquire mass through radiative corrections. Thus, we need a "hidden" sector 
where SUSY is broken spontaneously at tree level and then communicated 
to the MSSM by some "messengers" . 

There are two problems here: how is SUSY broken in the hidden sec- 
tor at the right scale and what are the messengers ? There are models in 
which a dynamical superpotential is generated by non-perturbative effects 
which breaks SUSY 0. The SUSY breaking scale is related to the Planck 
scale by dimensional transmutation. Two possibilities have been discussed 
in the literature for the messengers. One is gravity which couples to both 
the sectors 0. In a supergravity theory, there are non-renormalizable cou- 
plings between the two sectors which generate soft SUSY breaking operators 
in the MSSM once SUSY is broken in the "hidden" sector. In the absence 
of a flavor symmetry, this theory has to be fine tuned to give almost de- 
generate squarks and sleptons of the first two generations which is required 
by Flavor Changing Neutral Current (FCNC) phenomenology p, |^. The 
other messengers are the SM gauge interactions 0. In these models, the 
scalars of the first two generations are naturally degenerate since they have 
the same gauge quantum numbers. This is an attractive feature of these 
models, since the FCNC constraints are naturally avoided and no fine tuning 
between the masses of the first two generation scalars is required. If this 
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lack of fine tuning is a compelling argument in favour of these models, then 
it is important to investigate whether other sectors of these models are fine 
tuned. In fact, we will argue (and this is also discussed in |lT], |l2l) that 
the minimal model (to be defined in section 2) of low energy gauge mediated 
SUSY breaking requires a minimum 7% fine tuning to generate a correct 
vacuum {Z mass). Further, if a gauge-singlet is introduced to generate the 
"/i" and terms, then the minimal model of low energy gauge mediated 

SUSY breaking requires a minimum 1% fine tuning to correctly break the 
electroweak symmetry. These fine tunings makes it difficult to understand, 
within the context of these models, how SUSY is to offer some understanding 
of the origin of electroweak symmetry breaking and the scale of the Z and 
W gauge boson masses. 

Our paper is organized as follows. In section |^, we briefly review both the 
"messenger sector" in low energy gauge mediated SUSY breaking models that 
communicates SUSY breaking to the Standard Model and the pattern of the 
sfermion and gaugino masses that follow. Section |^ quantifles the flne tuning 
in the minimal model using the Barbieri-Giudice criterion 0. We show that 
a flne tuning of ~ 7% is required in the Higgs sector to obtain rriz- Section 
^ describes a toy model with split (5 + 5) messenger representations that 
improves the flne tuning. To maintain gauge coupling uniflcation, additional 
triplets are added to the low energy theory. They acquire a mass of 0(500) 
GeV by a coupling to a singlet. The flne tuning in this model is improved to 
~ 40%. The sparticle phenomenology of these models is also discussed. In 
section |, we discuss a version of the toy model where the above mentioned 
singlet generates the and terms. This is identical to the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) |]13| with a particular pattern for 
the soft SUSY breaking operators that follows from gauge mediated SUSY 
breaking and our solution to the flne tuning problem. We show that this 
model is tuned to ~ 20%, even if LEP does not discover SUSY/light Higgs. 
We also show that the NMSSM with one complete messenger (5 + 5) is flne 
tuned to ~ 1%. We discuss, in section |^, how it is possible to make our 
toy model compatible with a Grand Unifled Theory (GUT) [|14| based upon 
the gauge group SU{5) x SU{5). The doublet-triplet splitting mechanism of 
Barbieri, Dvali and Strumia [|1^] is used to split both the messenger represen- 
tations and the Higgs multiplets. In section 7, we present a model in which 
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all operators consistent with symmetries are present and demonstrate that 
the low energy theory is the model of section ^. In this model -R-parity {Rp) 
is the unbroken subgroup of a Z4 global discrete symmetry that is required to 
solve the doublet-triplet splitting problem. Our model has some metastable 
particles which might cause a cosmological problem. In the appendix, we 
give the expressions for the Barbieri-Giudice parameters (for the fine tuning) 
for the MSSM and the NMSSM. 



2 Messenger Sector 

In the models of low energy gauge mediated SUSY breaking |16[ (hence- 
forth called LEGM models), SUSY breaking occurs dynamically in a "hid- 
den" sector of the theory at a scale Adyn that is generated through dimen- 
sional transmutation. SUSY breaking is communicated to the Standard 
Model fields in two stages. First, a non-anomalous U{1) global symmetry of 
the hidden sector is weakly gauged. This U{l)x gauge interaction communi- 
cates SUSY breaking from the original SUSY breaking sector to a messenger 
sector at a scale Amess ~ c^xAdyn/ (4vr) as follows. The particle content in the 
messenger sector consists of fields 0+, 0_ charged under this U{l)x, a gauge 
singlet field S, and vector-like fields that carry Standard Model quantum 
numbers (henceforth called messenger quarks and leptons). In the minimal 
LEGM model, there is one set of vector-like fields, q, I, and q, I that together 
form a (5 + 5) of SU{5). This is a suffucient condition to maintain unification 
of the SM gauge couplings. The superpotential in the minimal model is 

Wmess = A^0+0_5 + hsS'' + X.Sqq + XiSll. (1) 

The scalar potential is 

V = J2\F,\' + ml\<P+\' + ml\<P_\'. (2) 



In the models of ||T0|, |l6l, the 0+, 0_ fields communicate (at two loops) with 



the hidden sector fields through the U{1) gauge interactions. Then, SUSY 
breaking in the original sector generates a negative value ~ — (axA^y^)^ / {4tt) 
for the mass parameters m^, of the 0+ and 0_ fields. This drives vevs of 
O (Afness) for the scalar components of both 0+ and 0_, and also for the scalar 



3 



and F-component of S if the couplings A5, gx and A</, satisfy the inequahties 
derived in |[Tl|, [1^.0 Generating a vev for both the scalar and F-component 



of S is crucial, since this generates a non-supersymmetric spectrum for the 
vector-like fields q and /. The spectrum of each vector- like messenger field 
consists of two complex scalars with masses ± B and two Weyl fermions 
with mass M where M = XS, B = XFs and A is the coupling of the vector- 
like fields to S. Since we do not want the SM to be broken at this stage, 
— B >0. In the second stage, the messenger fields are integrated out. 
As these messenger fields have SM gauge interactions, SM gauginos acquire 
masses at one loop and the sfermions and Higgs acquire soft scalar masses at 
two loops The gaugino masses at the scale at which the messenger fields 
are integrated out, A^ess ~ M are [|l^ 

M.^-(|==^WE^g(™)/. (3) 

The sum in equation |] is over messenger fields (m) with normalization 
Tr (T'^T^) = N^{m)6°'^ where the T's are the generators of the gauge group 
G in the representation R, fi{x) = 1 + 0(x), and Asusy = B/M = Fs/S = 
xAmess with X = B / M"^ . Henceforth, we will set Asusy ~ Amess- The 
exact one loop calculation []TB| of the gaugino mass shows that /i(a;) < 1.3 



for X < 1. The soft scalar masses at A„ 

= 2Al,sY E MC^f (^.) ( ^^^^^^ ) /2 ( ) , (4) 

where C'j^{si) is the Caismir of the representation of the scalar i in the gauge 




group G and f2{,x) = 1 + 0{x). The exact two loop calculation |]T8l which 
determines /2 shows that for x <0.8 (0.9), /2 differs from one by less than 
1%(5%). Henceforth we shall use fi{x) = 1 and f2{x) = 1. In the minimal 
LEGM model ^ 

47r 



^ This point in field space is a local minimum. There is a deeper minimum where SM 
is broken To avoid this problem, we can, for example, add another singlet to the 

messenger sector This does not change our conclusions about the fine tuning. 

^If all the dimensionless couplings in the superpotential are of 0(1), then x cannot be 
much smaller than 1. 
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m\K^ess) = 2A^,,, X (6) 

^ / «3(A mess mess 

)\ 3 /q;i(A 

mess 

^ V 47r / V 47r / 5 1^ An ) 

where Q = T^l + ^ and ai is the Sf/ (5) normahzed hypercharge couphng. 
Further, C3 = 4/3 and C2 = 3/4 for colored triplets and electroweak doublets 
respectively. 

The spectrum in the models is determined by only a few unknown param- 
eters. As equations 3 and 4 indicate, the SUSY breaking mass parameters 
for the Higgs, sfermions and gauginos are 

irig, m-g : tuh,, : me^,m^ ~ : as : ai. (7) 

The scale of Amess is chosen to be ~ 100 TeV so that the lightest of these par- 
ticles escapes detection. The phenomenology of the minimal LEGM model 



is discussed in detail in 19 . 



3 Fine Tuning in the Minimal LEGM 

A desirable feature of gauge mediated SUSY breaking is the natural sup- 
pression of FCNC processes since the scalars with the same gauge quantum 
numbers are degenerate 0. But, the minimal LEGM model introduces a 
fine tuning in the Higgs sector unless the messenger scale is low. This has 



been previously discussed in [10, 11| and quantified more recently in |12]. We 



outline the discussion in order to introduce some notation. 
The superpotential for the MSSM is 

W = flH^Hd + WYukan,a. (8) 

The scalar potential is 

V = ijllHuf + /isl^^r - ifJ'lHuHd + /i.c.)+D-terms + Vnoop, (9) 

where Vi-ioop is the one loop effective potential. The vev of Hu {Ha), denoted 
by Vu{vd), is responsible for giving mass to the up (down)-type quarks, /if = 
Tn'jj^ + A*^) f^2 — ^'hu A*3' ^'huI ^"fid SUSY breaking mass 

is often written as B/i. 
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terms for the Higgs fields. Q Extremizing this potential determines, with 
tan/3 = Vu/vd, 

= tan^/3-1 ' 

sin2/3 = 2 (11) 

where /i^ = fi] + 2dVi^ioop/ dvf . For large tan/3, m|/2 ^ + yU^)- This 

indicates that if |?tJ'//^| is large relative to m|, the yU^ term must cancel this 
large number to reproduce the correct value for m|. This introduces a fine 
tuning in the Higgs potential, that is naively of the order m|/(2|m|^^|). We 
shall show that this occurs in the minimal LEGM model. 

In the minimal LEGM model, a specification of the messenger particle 
content and the messenger scale A^ess fixes the sfermion and gaugino spec- 
trum at that scale. For example, the soft scalar masses for the Higgs fields are 
~ a2(Amess)Amess/(47r). Reuormalization Group (RG) evolution from Kmess 
to the electroweak scale reduces \m\j^ \ due to the large top quark Yukawa cou- 
pling, Af, and the squark soft masses. The one loop Reuormalization Group 
Equation (RGE) for m^j^ is (neglecting gaugino and the trilinear scalar term 
{HuQu'^) contributions ) 



dt 87r2 

which gives 



<At)+mUt)+mUt)), (12) 



(13) 

On the right-hand side of equation O the RG scaling of and m?c has 
been neglected. Since the logarithm \t\ ~ln(Amess/'"^t) is small, it is naively 
expected that m'j^^ will not be driven negative enough and will not trigger 
electroweak symmetry breaking. However since the squarks are ~ 500 GeV 
(1 TeV) for a messenger scale A^ess = 50 TeV (100 TeV), the radiative 
corrections from virtual top squarks are large since the squarks are heavy. 
A numerical solution of the one loop RGE (including gaugino and the tri- 
linear scalar term {H^Qu'^) contributions) determines —fn'n^ =(275 GeV)^ 

^The scale dependence of the parameters appearing in the potential is implicit. 
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((550 GeV)2) for A^^,, =50 TeV (100 TeV) and setting A* = 1. Therefore, 
m|/(2|m|^^|) ~0.06 (0.01), an indication of the fine tuning required. 

To reduce the fine tuning in the Higgs sector, it is necessary to reduce 
|m^^|; ideally so that m|^^ ^ — 0.5m|. The large value of |^|/^| at the weak 
scale is a consequence of the large hierarchy in the soft scalar masses at the 
messenger scale: < m^^ <^ ^'^ac- Models of sections and |^ attempt 
to reduce the ratio m^/m|^^ at the messenger scale and hence improve the 
fine tuning in the Higgs sector. 

The fine tuning may be quantified by applying one of the criteria of P, |]. 
The value O* of a physical observable O will depend on the fundamental 
parameters (Aj) of the theory. The fundamental parameters of the theory are 
to be distinguished from the free parameters of the theory which parameterize 
the solutions to 0{Xi) = O*. If the value O* is unusually sensitive to the 
underlying parameters (Aj) of the theory then a small change in Aj produces 
a large change in the value of O. The Barbieri-Giudice function 

quantifies this sensitivity 0. This particular value of O is fine tuned if the 
sensitivity to Aj is larger at O = O* than at other values of O [Q. If there 
are values of O for which the sensitivity to Aj is small, then it is probably 
sufficient to use c(0, Aj) as the measure of fine tuning. 

To determine c(m|, Aj), we performed the following. The sparticle spec- 
trum in the minimal LEGM model is determined by the four parameters 
^mess, /^35 and tan/3. The scale Amess fixes the boundary condition 
for the soft scalar masses, and an implicit dependence on tan/5 from A^, A;, 
and A^ arises in RG scaling^ from fi^c = Amess to the weak scale, that is 
chosen to be fij^Q = + \{jh1 + m^c). The extremization conditions of the 
scalar potential (equations |T0| and |11]) together with mz and leave two 
free parameters that we choose to be A^ess and tan (3 (see appendix for the 
expressions for these functions). 

A numerical analysis yields the value of c(m|,/i^) that is displayed in 
figure 1 in the (tan/5, A^ess) plane. We note that c(m|,/i^) is large through- 
out most of the parameter space, except for the region where tan /5 ^ 5 and 

®We allow for an arbitrary /i| at Kmess- 
^The RG scaling of At was neglected. 
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the messenger scale is low. A strong constraint on a lower limit for Amess 
is from the right-handed selectron mass. Contours mg^ = 75 GeV (~ the 
LEP limit from the run at ^ 170 GeV ^) and 85 GeV (~ the ultimate 
LEP2 limit [^]) are also plotted. The (approximate) limit on the neutralino 
masses from the LEP run at ^/s ~ 170 GeV, m^o + m^o = 160 GeV and the 
ultimate LEP2 limit, m^o +m^o ~ 180 GeV are also shown in figures a and c 
for sgn{fi) = —1 and figures b and d for sgn{fi) = +1. The constraints from 
the present and the ultimate LEP2 limits on the chargino mass are weaker 
than or comparable to those from the selectron and the neutralino masses 
and are therefore not shown. If rriz were much larger, then c ~ 1. For ex- 
ample, with mz = 275 GeV (550 GeV) and Amess= 50 (100) TeV, c(m|;//2) 
varies between 1 and 5 for 1.4 ^ tan/5 ^ 2, and is ~ 1 for tan/3 > 2. This 
suggests that the interpretation that a large value for c(m|; /x^) implies that 
TUz is fine tuned is probably correct. 

From figure 1 we conclude that in the minimal LEGM model a fine tuning 
of approximately 7% in the Higgs potential is required to produce the correct 
value for mz- Further, for this fine tuning the parameters of the model are 
restricted to the region tan/? ^ 5 and A^ess ~ 45 TeV, corresponding to 
mg^ ^ 85 GeV. We have also checked that adding more complete (5 + 5)'s 
does not reduce the fine tuning. 



4 A Toy Model to Reduce Fine Tuning 

4.1 Model 

In this section the particle content and couplings in the messenger sector that 
are suffucient to reduce is discussed. The aim is to reduce m'^/m'jj^ 

at the scale A^g^^. 

The idea is to increase the number of messenger leptons {SU (2) doublets) 
relative to the number of messenger quarks {SU{3) triplets). This reduces 
both m'^/m'jj and m%/ml at the scale Amess (see equation Hf). This leads 
to a smaller value of in the RG scaling (see equation 0) and the scale 

Amess can be lowered since mg^ is larger. For example, with three doublets 
and one triplet at a scale A^ess = 30 TeV, so that mg^ ~ 85 GeV, we find 
|m^^(mQ)| ^ (lOOGeV)^ for Xt = 1. This may be achieved by the following 
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superpotential in the messenger sector 

W = X,,Sqiqi + Xi,Shh + Xi,Sl2l2 + KShl3 + ^XsS' 

+X^S^^4>+ + ^XnN^ + Xg,Nq2q2 + Xg^Nq^q^, (15) 

where iV is a gauge singlet. The two pairs of triplets ^2,^2 and 53,^3 are 
required at low energies to maintain gauge coupling unification. In this model 
the additional leptons /2,^2 and /3,r3 couple to the singlet S, whereas the 
additional quarks couple to a different singlet that does not couple to the 
messenger fields 0+, 0_. This can be enforced by discrete symmetries (we 
discuss such a model in section |^). Further, we assume the discrete charges 
also forbid any couplings between N and S at the renormalizable level (this 
is true of the model in section 0) so that SUSY breaking is communicated 
first to S and to N only at a higher loop level. 

4.2 Mass Spectrum 

Before quantifying the fine tuning in this model, the mass spectrum of the ad- 
ditional states is briefiy discussed. While these fields form complete represen- 
tations of SU{5), they are not degenerate in mass. The vev and F-component 
of the singlet S gives a mass Amess to the messenger lepton multiplets if the 
F-term splitting between the scalars is neglected. As the squarks in -|- (ji 
(i=2,3) do not couple to S, they acquire a soft scalar mass from the same 
two loop diagrams that are responsible for the masses of the MSSM squarks, 
yielding rrig ~ a3{Amess) Asusy /{VQt^)- The fermions in g + g also acquire 
mass at this scale since, if either or ~ 0(1), a negative value for 
m% (the soft scalar mass squared of A^) is generated from the XgNqq cou- 
pling at one loop and thus a vev for A^ ~ rrig is generated. The result is 
mi/rrig^ V67r/a3(A 

mess )(A mess /Asusy) ~ 85. 
The mass splitting in the extra fields introduces a threshold correction 
to sin^ 6*14^ if it is assumed that the gauge couplings unify at some high scale 
Mgut ~10^^ GeV. We estimate that the splitting shifts the prediction for 
siv? Ow by an amount ~ — 7x 10""^ ln(m;/mg)n, where n is the number of 
split (5 + 5).|3 In this case n =2 and mi/nig ~ 85, so 5sin^ 6'vi/ ~ — 6 x 
^"The complete (5 + 5), i.e., hji and qi,qi, that couples to S is also split because 
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10 ^. If a^^Mz) and aem{Mz) are used as input, then using the two loop 
RG equations sm^OwiMS) = 0.233 ± 0(10-^) is predicted in a minimal 



SUSY-GUT [221. The error is a combination of weak scale SUSY and GUT 



threshold corrections!^. The central value of the theoretical prediction is 
a few percent higher than the measured value of sin^ ^^(MS*) = 0.231 ± 
0.0003|^. The split extra fields shift the prediction of sin^6'iy to ~ 0.227 ± 



O(10~ ) which is a few percent lower than the experimental value. In sections 
we show that this spectrum is derivable from a SU{5) x SU{5) GUT 
in which the GUT threshold corrections to sin^ could be ~ 0(10^'^) — 
O(10~^) 12^. It is possible that the combination of these GUT threshold 
corrections and the split extra field threshold corrections make the prediction 
of sin^ 9w niore consistent with the observed value. 

4.3 Fine Tuning 

To quantify the fine tuning in these class of models the analysis of section 3 
is applied. In our RG analysis the RG scaling of At, the effect of the extra 
vector-like triplets on the RG scaling of the gauge couplings, and weak scale 
SUSY threshold corrections were neglected. We have checked a postiori that 
this approximation is consistent. As in section ^, the two free parameters are 
chosen to be Amess and tan/5. Contours of constant c(m|,/i^) are presented 
in figure 2. We show contours of m^o + m^o = 160 GeV, and = 75 
GeV in figure 2 a for sgn{fi) = —1 and 2b for sgn{fi) = +1. These are 
roughly the present limits from LEP (including the run at ^/s ~ 170 GeV 
2^). The (approximate) ultimate LEP2 reaches m^o +m^o = 180 GeV, 



and mg^ = 85 GeV are shown in figure 2c for sgn{fi) = — 1 and figure 2d 
for sgn{fi) = +1. Since /i^(~ (100 GeV)^) is much smaller in these models 
than in the minimal LEGM model, the neutralinos (x? and xt) lighter 
so that the neutralino masses provide a stronger constraint on Amess than 
does the slepton mass limit. The chargino constraints are comparable to the 
neutralino constraints and are thus not shown. It is clear that there are areas 
of parameter space in which the fine tuning is improved to ~ 40% (see figure 



A( 7^ Xq at the messenger scale due to RG scaling from Mqut to Amess- This splitting is 
small and neglected. 
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While this model improves the fine tuning required of the jj, parameter, it 
would be unsatisfactory if further fine tunings were required in other sectors 
of the model, for example, the sensitivity of m| to Amess and and the 
sensitivity of rrit to /i^, Amess and Xt. We have checked that all these are 
less than or comparable to c{m%; /x^). We now discuss the other fine tunings 
in detail. 

For large tan/?, the sensitivity of m| to c(m|; cx 1/ tan^/3, and 
is therefore smaller than c(m|;/i^). Our numerical analysis shows that 
c(m|; nl) ^ c(m|; /x^) for all tan/3. 

In the one loop approximation m|^^ and m|^^ at the weak scale are pro- 
portional to A^gj,^ since all the soft masses scale with A^ess and there is 
only a weak logarithmic dependence on Amess through the gauge couplings. 
We have checked numerically that {A^gg/m\j^){dm%^/ dA^^^ ~ 1. Then, 
c(^l;^meJ ~ c(m|;m|^J +c(m|;m|^J. We find that c(m|;A^gJ 
c{m?2:] over most of the parameter space. 

In the one loop approximation, m\^(t) is 

m|Ji) ^ m|jO) + (m|3(0) +m?c(0) +m^jO))(e-i^* - 1). (16) 



Then, using t fti hi[Amess/^Q^) ~ ln(\/67r/Q;3) 4.5 and At ^ 1, c(m|; At) is 
(see appendix) 

, . , 4 drriir it) 

This result measures the sensitivity of m| to the value of At at the electroweak 
scale. While this sensitivity is large, it does not reflect the fact that Xt{Mpi) 
is the fundamental parameter of the theory, rather than Xt{M^eak)- We 
flnd by both numerical and analytic computations that, for this model with 
three (5 + 5)'s in addition to the MSSM particle content, SXt{Myjeak) ~ 
0.1 X SXt{Mpi), and therefore 

2 

in ~ 

c(m|;A,(M„))«5j^;j^. (18) 

For a scale of Amess = 50 TeV (m^ ^ 600 GeV), c(m|; At(Mp;)) is comparable 
to c(m|; /x^) which is ~ 4 to 5. At a lower messenger scale, Amess ~ 35 
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rrif „ rrizc 


687 


616 


612 


319 125 





rriuc 


656 


546 



Table 1: Soft scalar masses in GeV for messenger particle content of three 
{I + lys and one q + q and a scale A^ess — 50 TeV. 

TeV, corresponding to squark masses of 450 GeV, the sensitivity of m| 
to Xt{Mpi) is 2.8. This is comparable to c(m|;//^) evaluated at the same 
scale. 

We now discuss the sensitivity of nit to the fundamental parameters. 
Since, — sin^ /3Af , we get 

c{mt; Xi) = 5x,Xi + ^c{mz; Xi) + ^.^^ (19) 

Numerically we find that the last term in c(mt; Xi) is small compared to 
c(m|; Xi) and thus over most of parameter space c{mt; Xi) ~ ^c{m\\ Xi). As 
before, the sensitivity of uit to the value of at the GUT/Planck scale is 
much smaller than the sensitivity to the value of Xt at the weak scale. 

4.4 Sparticle Spectrum 

The sparticle spectrum is now briefly discussed to highlight deviations from 
the mass relations predicted in the minimal LEGM model. For example, 
with three doublets and one triplet at a scale of A = 50 TeV, the soft scalar 
masses (in GeV) at a renormalization scale = "^t + |("^|3+"^Ig) ~ (630 
GeV)^, for At = 1, are shown in table 1. 

Two observations that are generic to this type of model are: (i) By con- 
struction, the spread in the soft scalar masses is less than in the minimal 
LEGM model, (ii) The gaugino masses do not satisfy the one-loop SUSY- 
GUT relation Mj/aj = constant. In this case, for example, M^/a^ : M2/a2 ~ 
1:3 and M^/a^ : Mi/ai ^ 5:11 to one-loop. 
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We have also found that for tan /5 3, the Next Lightest Supersymmetric 
Particle (NLSP) is one of the neutralinos, whereas for tan (3 ^ 3, the NLSP 
is the right-handed stau. Further, for these small values of tan (3, the three 
right-handed sleptons are degenerate within 200 MeV. 



5 NMSSM 

In section ^, the fi term and the SUSY breaking mass fi^ were put in by 
hand. There it was found that these parameters had to be fine tuned in 
order to correctly reproduce the observed Z mass. The extent to which this 
is a "problem" may only be evaluated within a specific model that generates 
both the fj, and /i| terms. 

For this reason, in this section a possible way to generate both the /i term 
and term in a manner that requires a minimal modification to the model 
of either section ^ or section ^ is discussed. The easiest way to generate these 
mass terms is to introduce a singlet N and add the interaction NHuH^ to 
the superpotential (the NMSSM) p|. The vev of the scalar component of 



generates /i and the vev of the F-component of N generates 

We note that for the "toy model" solution to the fine tuning problem 
(section the introduction of the singlet occurs at no additional cost. Recall 
that in that model it was necessary to introduce a singlet A^, distinct from 
S, such that the vev of gives mass to the extra light vector-like triplets, 



qi,qi {i = 2,3) (see equation |T5|) . Further, discrete symmetries (see section 
1^) are imposed to isolate from SUSY breaking in the messenger sector. 
This last requirement is necessary to solve the fine tuning problem: if both 
the scalar and F-component of A^ acquired a vev at the same scale as S, then 
the extra triplets that couple to A^ would also act as messenger fields. In this 
case the messenger fields would form complete (5 + 5)'s and the fine tuning 
problem would be reintroduced. With A^ isolated from the messenger sector 
at tree level, a vev for A^ at the electroweak scale is naturally generated, as 
discussed in section ||. 

We also comment on the necessity and origin of these extra triplets. Re- 
call that in the toy model of section ^ these triplets were required to maintain 
the SUSY-GUT prediction for sin^ 6w- Further, we shall also see that they 
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are required in order to generate a large enough —mj^ (the soft scalar mass 
squared of the singlet N). Finally, in the GUT model of section the light- 
ness of these triplets (as compared to the missing doublets) is the consequence 
of a doublet-triplet splitting mechanism. 

The superpotential in the electroweak symmetry breaking sector is 

W = + \Nqq - XhNH^H,, (20) 

which is similar to an NMSSM except for the coupling of to the triplets. 
The superpotential in the messenger sector is given by equation |15|. 
The scalar potential is |^ 

V = Y.\^i\^ + ^li\N\^ + ^HjHu\^ + mjjjHd\^ + 'D-teTms 

i 

-{AhNH^H^ + h.c.) + V.^ioop. (21) 

The extremization conditions for the vevs of the real components of A^, 
and Hd, denoted by vn, Vu and Vd respectively (with v = \fv'^ + Vd ~ 250 
GeV), are 

VN{rh% + Xh— + X%v% - XHXNVuVd) i=^i?f«fd = 0, (22) 



with 



~ tan2/3-l ' ^^^^ 



sin 2/3 = 2 (24) 



1^' = (25) 

11 1 

(J'l = --XjjVuVd+ -XhXnv% + Ah-j^vn, (26) 

= + i = {u,d,N). (27) 



^^In models of gauge mediated SUSY breaking, Ah=0 at tree level and a non-zero value 
of Ah is generated at one loop. The trilinear scalar term AnN^ is generated at two loops 
and is neglected. 
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We now comment on the expected size of the Yukawa couphngs Xg, Xn 
and Xh- We must use the RGE's to evolve these couphngs from their val- 
ues at Mgut or Mpi to the weak scale. The quarks and the Higgs doublets 
receive wavefunction renormalization from S't/(3) and SU{2) gauge interac- 
tions respectively, whereas the singlet does not receive any wavefunction 
renormalization from gauge interactions at one loop. So, the couplings at 
the weak scale are in the order: ~ 0(1) > Xh > Xn if they all are 0(1) 
at the GUT/Planck scale. 

We remark that without the Nqq coupling, it is difficult to drive a vev 
for as we now show below. The one loop RGE for is 

(28) 

Since A^ is a gauge- singlet, = at A^ess- Further, if Xg = 0, an estimate 
for m% at the weak scale is then 

ml - -|f (-L(O) + <(0)) In , (29) 

i.e., Xh drives negative. The extremization condition for vn, equation 



|, and using equations ^ and (neglecting Ah) shows that 



m 



^ + Al,^ - XI f ^ - ^(mijO) + m^^(O)) In ( ^]] (30) 



87r2^ ' " V mn 



d 

has to be negative for A^ to acquire a vev. This implies that m|^^ and m|^^ at 
^mess have to be greater than ~ (350 GeV)^ which implies that a fine tuning 
of a few percent is required in the electroweak symmetry breaking sector. 
With Xg ~ 0(1), however, there is an additional negative contribution to 
mj^ given approximately by 

-SK(0)+^(0))lnf^). (31) 



This contribution dominates the one in equation ^ since Xg > Xh and the 
squarks q, q have soft masses larger than the Higgs. Thus, with Xg ^ 0, 



m 



N 



+ A^f2/2 is naturally negative. 



Fixing mz and rrit, we have the following parameters : Amess, Xg, Xh, Xn, 
tan/?, and v^. Three of the parameters are fixed by the three extremization 
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conditions, leaving three free parameters that for convienence are chosen to 
be A^ess, tan/3 >0, and Xh- The signs of the vevs are fixed to be positive 
by requiring a stable vacuum and no spontaneous CP violation. The three 
extremization equations determine the following relations 

Aiv = j^{fil + ^Xjjsm2f3v^-^AHVN), (32) 

Vn = V2^, (33) 

1 1 1 

ml = \N\H7^sm2Pv^-X%v%~-\lv^ + —^AHsm2P—, (34) 

2 2 2v2 Vn 



where 



^ = "2"^^+ 1-tan^/? ' 
2fil = sm2P{2fi' + mjj^+mlJ. (36) 



m 



N- 



The superpotential term NHuHd couples the RGB's for mj^^, ttih^ and 
Thus the values of these masses at the electroweak scale are, in general, 
complicated functions of the Yukawa parameters At, A^, Aat and Xq. In our 
case, two of these Yukawa parameters (A^ and A^r) are determined by the 
extremization equations and a closed form expression for the derived quan- 
tities cannot be found. To simplify the analysis, we neglect the dependence 
of Trijj^ and m|^^ on Xfj induced in RG scaling from A^ess to the weak scale. 
Then and rrifj^ depend only on Amess and tan/5 and thus closed form 
solutions for Aat, vn and can be obtained using the above equations. 
Once rh% at the weak scale is obtained, the value of Xq is obtained by using 
an approximate analytic solution. An exact numerical solution of the RGE's 
then shows that the above approximation is consistent. 



5.1 Fine Tuning and Phenomenology 

The fine tuning functions we consider below are c(0; A^), c(0; Aat), c(0; Xt), 
c(0; Xq) and c(0; Amess) where O is either m| or m^. The expressions for the 
fine tuning functions and other details are given in the appendix. In our RG 
analysis the approximations discussed in subsection ^]3| and above were used 
and found to be consistent. Fine tuning contours of c(m|; Xh) are displayed 
in figures 3 a and 3 b for Xh = 0.1 and figures 3 c and 3 d for Xh = 0.5. We 
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have found by numerical computations that the other fine tuning functions 
are either smaller or comparable to c(m|; Aj^). 

We now discuss the existing phenomenological constraints on our model 
and also the ultimate constraints if LEP2 does not discover SUSY/light 
Higgs(/i). These are shown in figures 3 a, 3 c and figures 3 b, 3 d respec- 
tively. We consider the processes e"'"e~ — > Zh^ e~^e~— {h + pseudoscalar) , 
e~^e~-^ X^X~y e"'"e~— > X1X2) e^e~—^ (^r(^*r observable at LEP. Since 
this model also has a light pseudoscalar, we also consider upsilon decays 
(7 + pseudoscalar) . We find that the model is phenomenologically vi- 
able and requires a ~ 20% tuning even if no new particles are discovered at 
LEP2. 

We begin with the constraints on the scalar and pseudoscalar spectra of 
this model. There are three neutral scalars, two neutral pseudoscalars and 
one complex charged scalar. We first consider the mass spectrum of the 
pseudoscalars. At the boundary scale Amess, SUSY is softly broken in the 
visible sector only by the soft scalar masses and the gaugino masses. Further, 



the superpotential of equation 20 has an /^-symmetry. Therefore, at the tree 



level, i.e., with Ah =0, the scalar potential of the visible sector (equation 
21| ) has a global symmetry. This symmetry is spontaneously broken by the 
vevs of A^'^, H^, and H^^ (the superscript R denotes the real component of 
fields), so that one physical pseudoscalar is massless at tree level. It is 

(vnN^ + V sin 2(3 cos (3HI + v sin 2(3 sin (3Hpj , (37) 



sin^ 2p 



where the superscripts / denote the imaginary components of the fields. The 
second pseudoscalar, 

A ^ -±N' + + (38) 

vtv fsmp fcosp 

acquires a mass 

m\ = -XHXNvlfitan (3 + cot P) + XhXnv"^ sin 2/3 (39) 



In computing these functions the weak scale value of the couplings Xn and Xh 
has been used. But since Xn and Xh do not have a fixed point behavior, we have 
found that Xh{Mgut) / XHimz) dXH{mz) / dXHiMcur) ~ 1 so that, for example, 
c(to|; Xh{Mgut)) ~ c{m\\ Xnimz)). 
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through the I-FatP term in the scalar potential. 

The pseudoscalar a acquires a mass once an An-term is generated, at 
one loop, through interactions with the gauginos. Including only the wino 
contribution in the one loop RGE, Ah is given by 



A 



H 



6- . 
20 \h 



-M2XH In 
M2 



Mo 



GeV, 



(40) 



.280GeVy 

where M2 is the wino mass at the weak scale. Neglecting the mass mixing 
between the two pseudoscalars, the mass of the pseudo-Nambu-Goldstone 
boson is computed to be 

9 



AvNVuVd/ {vj^ + V sin 2/3) 




Mo 



( 



sin 2/5 



250GeV 



sm' 



(250GeV 



(GeV)2(41) 



280GeV 

If the mass of a is less than 7.2 GeV, it could be detected in the decay 



T 
T 



/i 



a + 7[|3 



Comparing the ratio of decay width for T — a + 7 to 



the limit 



sin 2/5 tan /3 



25ufev)^ + --^2/5 



< 0.43 



(42) 



is found. 

Further constraints on the spectra are obtained from collider searches. 
The non-detection oi Z ^ scalar + a at LEP implies that the combined 
mass of the lightest Higgs scalar and a must exceed ~ 92 GeV. Also, the 
process e^e~ ^Zh may be observable at LEP2. For \h = 0.1, the constraint 
rrifi + rria ^ 92 GeV is stronger than ^ 70 GeV which is the limit from 
LEP at ^/s ~ 170 GeV [20|. The contour of rrih + rria = 92 GeV is shown in 
figure 3 a. In figure 3 b, we show the contour of rrih = 92 GeV (~ the ultimate 
LEP2 reach ||26|). For Xh = 0.5, we find that the constraint rrih ^ 70 GeV 
is stronger than rrih + ^ 92 GeV and restricts tan/5 ^ 5 independent of 
Amess- The contour rrih = 92 GeV is shown in figure 3 d. We note that the 
allowed parameter space is not significantly constrained. We find that these 
limits make the constraint of equation redundant. The left-right mixing 



18 



between the two top squarks was neglected in computing the top squark 
radiative corrections to the Higgs masses. 

The pseudo-Nambu-Goldstone boson a might be produced along with the 
lightest scalar h at LEP. The (tree- level) cross section in units of i? = 87/s 
nb is 



a(e+e- ^ /i a) ^ 0.15- i; f 1, , (43) 

[s — miY \ s s 



where g\/ cos6'iy is the Z{adh — hda) coupling, and 

f (x, y, z) = \J {x — y — zY — Ayz. \i h = cnN^ + CuH^ + CdH^, then 

A = sin 2/3 "f<^-^'°fe . (44) 

We have numerically checked the parameter space allowed by rrih ^ 70 GeV 
and \h <0.5 and have found the production cross section for ha to be less 
than both the current limit set by DELPHI [E^ and a (possible) exclusion 



limit of 30 fb at y/s ^ 192 GeV. The production cross-section for hA is 
larger than for ha and A is therefore in principle easier to detect. However, for 
the parameter space allowed by rrih ^ 70 GeV, numerical calculations show 
that rriA ^ 125 GeV, so that this channel is not kinematically accessible. 
The charged Higgs mass is 

ml± =mlr + m^^ + m^^ + 2fi^ (45) 

which is greater than about 200 GeV in this model since m|^^ ^ (200GeV)2 
for Amess ^ 35 TeV and as ~ — m|^^. 

The neutralinos and charginos may be observable at LEP2 at ^/s ~ 192 
GeV if ^ 95 GeV and m^o +m^o ^ 180 GeV. These two constraints are 

A, Xl A2 

comparable, and thus only one of these is displayed in figures 3 b and 3 d, 
for Ah = 0.1 and Xh = 0.5 repectively. Also, contours of m^o + m^o = 160 
GeV (~ the LEP kinematic limit at y/s ~ 170 GeV) are shown in figures 3 
a and 3 c. Contours of 85 GeV (~ the ultimate LEP2 limit) and 75 GeV 
(~ the LEP limit from y/s ~ 170 GeV) for the right-handed selectron mass 
further constrain the parameter space. 

The results presented in all the figures are for a central value of mt=175 
GeV. We have varied the top quark mass by 10 GeV about the central value 
of 1711= 175 GeV and have found that both the fine tuning measures and the 
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LEP2 constraints (the Higgs mass and the neutrahno masses) vary by ^ 30 
%, but the quahtative features are unchanged. 

We see from figure 3 that there is parameter space allowed by the present 
limits in which the tuning is ~ 30 %. Even if no new particles are discovered 
at LEP2, the tuning required for some region is ~ 20%. 

It is also interesting to compare the fine tuning measures with those found 
in the minimal LEGM model (one messenger (5 + 5)) with an extra singlet 

to generate the /i and /ig terms.|^ In figure 4 the fine tuning contours 
for c(m|; Ai^) are presented for Xh=0.1. Contours of = 75 GeV and 
m^o + m^o = 160 GeV are also shown in figure 4 a. For Xh = 0.1, the 
constraint m/j + nia ^ 92 GeV is stronger than the limit rrih ^ 70 GeV and 
is shown in the figure 4 a. In figure 4 b, we show the (approximate) ultimate 
LEP2 limits, i.e., nih = 92 GeV, m^o+m^o = 180 GeV and mg^ = 85 GeV. Of 
these constraints, the bound on the lightest Higgs mass (either mh + nia ^ 92 
GeV or rrih ^ 92 GeV) provides a strong lower limit on the messenger scale. 
We see that in the parameter space allowed by present limits the fine tuning 
is ^ 2% and if LEP2 does not discover new particles, the fine tuning will be 
^ 1%. The coupling Xh is constrained to be not significantly larger than 0.1 
if the constraint rrih + "^a ^ 92 GeV (or rrih ^ 92 GeV) is imposed and if the 
fine tuning is required to be no worse than 1%. 

6 Models Derived from a GUT 

In this section, we discuss how the toy model of section | could be derived 

from a GUT model. 

In the toy model of section ^, the singlets and 5* do not separately 

couple to complete SU{5) representations (see equation |I3p. If the extra 

fields introduced to solve the fine tuning problem were originally part of 

(5 + 5) multiplets, then the missing triplets (missing doublets) necessarily 

couple to the singlet S{N). The triplets must be heavy in order to suppress 

their contribution to the soft SUSY breaking mass parameters. If we assume 

the only other mass scale is Mqut, they must acquire a mass at Mqut- This 

^•^We assume that the model contains some mechanism to generate —m?j^ ^ (lOOGeV)^ — 
(200GeV)^; for example, the singlet is coupled to an extra (5 + 5). 



20 



is just the usual problem of splitting a (5 + 5) |]T4|. For example, if the 
superpotential in the messenger sector contains four (5 + 5)'s, 



W = XiSliihi + A255,25z2 + A355z35/3 + KSl^bg, (46) 

then the SU{?>) triplets in the (5; + 5;)'s and the SU{2) doublet in (5^ + 5g) 
must be heavy at Mqut so that in the low energy theory there are three 
doublets and one triplet coupling to S. This problem can be solved using 
the method of Barbieri, Dvali and Strumia [|l5l that solves the usual Higgs 
doublet-triplet splitting problem. The mechanism in this model is attractive 
since it is possible to make either the doublets or triplets of a quintet heavy 
at the GUT scale. We next describe their model. 

The gauge group is SU{b) x S'f/(5)', with the particle content S(24, 1), 
S'(l, 24), $(5, 5) and4>(5, 5) and the superpotential can be written as 

W = $;^,(M$5^;5^ + AS;^5°; + A'S';,'5;^)$f + 
+^MsTr(S2) + ^Ms'Tr(S'2) + 

^AsTrS^^ + ^As'IVE'3. (47) 

A supersymmetric minimum of the scalar potential satisfies the F - flatness 
conditions 

= F^ = {M^5p1 + \^5^; + \'T.'^:51)^^;, 

= Fs = lMsS^ + i(A<l>^'-A^5fTr(<l<|.))+AE(S2-^TrS2), 
= = ^Ms^sJ + i (A'<l^,<|.f - A'^52Tr($<|.)) + AsKS^^ 

(48) 

With the ansatz Q 

S = diag(2, 2, 2, -3, -3) , S' = v^, diag(2, 2, 2, -3, -3), (49) 
the = condition is 

diag[M3, M3, M3, M2, M2] ■ diagh, v^, v^, V2, V2] = 0, (50) 

^''The two possible solutions to the F-flatness conditions are — vs diag(2, 2, 2, —3, —3) 
and S = diag(l, 1, 1, 1, —4). 
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where M3 = M$ + 2Afs + 2A'fE' and M2 = M$ — 3Afs — 3A'fE' and the second 
matrix is the vev of $. To satisfy this condition, there is a discrete choice 
for the pattern of vev of $ : i) f 3 7^ and M3 = or ii) f 2 7^ and M2 = 0. 
Substituting either i) or ii) in the F^, and F-^i conditions then determines ^3 
(or ^2). With two sets of fields, $1, <l>i with f 3 7^ and $2, ^2 with V2 ^ , 
we have the following pattern of symmetry breaking 

SU{5) X SU{5y '''4^' {SU{3) X SU{2) x U{1)) x {SU{3) x SU{2) x f/(l))' 

S'M (the diagonal subgroup). (51) 

If the scales of the two stages of symmetry breaking are about equal, i.e. 
fs, v-£', ~ f 3 , f 2 ~ Mgut, then the SM gauge couplings unify at the scale 
Mgut. □ 

The particular structure of the vevs of $1 and $2 can be used to split 
representations as follows. 

Consider the Higgs doublet-triplet splitting problem. With the particle 
content 5/i(5, 1), 5^(5, 1) and X(l,5), X(l,5) and the superpotential 

W = 5;,„X"'<I?„, + 5^X„,$i°', (52) 

the SU{3) triplets in 5h, 5h and X, X acquire a mass of order Mgut whereas 
the doublets in 5h, 5h and X, X are massless. We want only one pair of 
doublets in the low energy theory (in addition to the usual matter fields). 
The doublets in X, X can be made heavy by a bare mass term Mgut^^ ■ 
Then the doublets in 5h, 5h are the standard Higgs doublets. But if all terms 
consistent with symmetries are allowed in the superpotential, then allowing 
Mgut^i^i, MgutXX , 5hX^i and 5hX^i implies that a bare mass term 
for 5h5h is allowed. Of course, we can by hand put in a /x term fj,5h5h of 
the order of the weak scale as in section ^. However, it is theoretically more 
desirable to relate all electroweak mass scales to the original SUSY breaking 
scale. So, we would hke to relate the /i term to the SUSY breaking scale. 
We showed in section ^ that the NMSSM is phenomenologically viable and 
"un-fine tuned" in these models. 

^^See and ||2^ for models which give this structure of vevs for the <I> fields without 
using the adjoints. 
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The vev structure of $2, ^'2 can be used to make the doublets in a 5 + 5 
heavy. Again, we get two pairs of hght triplets and one of these pairs can be 
given a mass at the GUT scale. 

We can use this mechanism of making either doublets or triplets in a 
(5 + 5) heavy to show how the model of section ^ is derivable from a GUT. 
The model with three messenger doublets and one triplet is obtained from a 
GUT with the following superpotential 

W = S55 + S5i5i + SX1X1 + 
5iXi^i + 5/X/$i + 

McurXhXh + 5hXh^^ + 5hXh<l>i + fi5h5h 

+Ar3 + N5q5g + NXqXq. (53) 

Here, some of the "extra" triplets and doublets resulting from splitting (5 + 
5)'s are massless at the GUT scale. For example, the "extra" light doublets 
are used as the additional messenger leptons. After inserting the vevs and 
integrating out the heavy states, this corresponds to the superpotential in 



equation |15| with the transcription: 



5,5 gi,gi + /i,/i 

5/, 5z I2, k 

XuXi — > Z3,Z3 

5q — > ^2 

Xq,Xq q3,q3. (54) 

We conclude this section with a remark about light singlets in SUSY- 
GUT's with low energy gauge mediated SUSY breaking.[^ In a SUSY GUT 
with a singlet coupled to the Higgs multiplets, there is a potential problem 
of destabilising the m^eak/MGUT hierarchy, if the singlet is light and if the 
Higgs triplets have a SUSY invariant mass of 0{Mgut) @- In the LEGM 
models, a B-type mass for the Higgs triplets and doublets is generated at 
one loop with gauginos and Higgsinos in the loop, and with SUSY breaking 
coming from the gaugino mass. Since SUSY breaking (the gaugino mass and 
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the soft scalar masses) becomes soft above the messenger scale, A„iess ~ 100 
TeV, the B-type mass term generated for the Higgs triplets is suppressed, i.e., 
it is 0((a/47r)M2ALss/^Gi/T)- Similarly the soft mass squared for the Higgs 
triplets are 0(m^g^^A^g^^/M^^j.). Since the triplets couple to the singlet N, 
the soft scalar mass and B-term generates at one loop a linear term for 
the scalar and F-component of respectively. These tadpoles are harmless 
since the SUSY breaking masses for the triplets are so small. This is to be 
contrasted with supergravity theories, where the 5-term~ 0{mweakMGUT) 
and the soft mass ~ 0{mu,eak) for the triplet Higgs generate a mass for the 
Higgs doublet that is at least ~ 0{y/mweakMGUT / (47r)). 

7 One complete Model 

The model is based on the gauge group Gioc = SU{5) x SU{5y and the 
global symmetry group Ggio = Z3 x Z'^ x Z4. The global symmetry acts 
universally on the three generations of the SM. The particle content and 
their Gioc x Ggio quantum numbers are given in table 2. The most general 
renormalizable superpotential that is consistent with these symmetries is 



W ^Wi + W2 + W3 + W4 + W5 + W6 + W7, 



(55) 



where. 




(56) 
(57) 
(58) 
(59) 

(60) 



W2 = M.XiX, 

W3 = Xi5h^iXh + Xi5h^iXh + X25i^iXi + X25i^iXi, 

= A35,$2^, + A35,$2^„ 





(61) 
(62) 



= Ag5a0,-5;, + AglOilO.-S,,. 



24 



The origin of each of the W^j's appearing in the superpotential is easy to 
understand. In computing the F=0 equations at the GUT scale, the only 
non-trivial contributions come from fields appearing in Wi , since all other WiS 
are bilinear in fields that do not acquire vevs at the GUT scale. The function 
of Wi is to generate the vevs S, S' ~ diag [2, 2, 2, -3, -3], = $2 ~ diag 
[0, 0, 0, 1, 1] and 5>f = $1 ~ diag [1, 1, 1, 0, 0]. These vevs are necessary to 
break Gioc -^SU{3)c x SU(2) x f/(l)y (this was explained in section H). The 
role of and W^^ is to generate the necessary splitting within the many 
(5 + 5)'s of Gioc that is necessary to solve the usual doublet-triplet splitting 
problem, as well as to solve the fine tuning problem that is discussed in 
sections |^,^ and ^ The messenger sector is given by W^. It will shortly 
be demonstrated that at low energies this sector contains three vector-like 
doublets and one vector-like triplet. The couplings in and Wj at low 
energies contain the electroweak symmetry breaking sector of the NMSSM, 
the Yukawa couplings of the SM fields, and the two light vector-like triplets 
necessary to maintain the few percent prediction for sin^ 6w as well as to 
generate a vev for A^. 

We now show that the low energy theory of this model is the model that 
is discussed in section |[ 

Inserting the vevs for $1 and $1 into W3, the following mass matrix for 
the colored triplet chiral multiplets is obtained: 



(5/i, 5;, Xi) 



( Aif$i 
Aif$i 




V 
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(63) 



and all other masses are zero. There are a total of four vector-like colored 
triplet fields that are massive at Mqut- These are the triplet components of 
(5/i,X/j), (5/i,X/j), (5/,Xi) and {Xi,Th), where Th is that linear combination 
of triplets in 5/ and X that marries the triplet component of X/. The orthog- 
onal combination to Th, Ti, is massless at this scale. The massless triplets 
at Mgut are (5^, 5g), {Xq,Xq) and (X,Ti), for a total of three vector-like 
triplets. By inspection, the only light triplets that couple to S* at a renor- 
malizable level are 5g and 5,, which was desirable in order to solve the fine 
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tuning problem. Further, since X contains a component of T^, the couphngs 
of the other hght triplets to the singlets N and A^' are 



Weff = XioN'XXg + XhN'X.Tl + X.NTlX, 



(64) 



where Xg = AgCosa', An = An cos a' and a' is the mixing angle between the 
triplets in 5/ and X, i.e., = cosa'X — sina'5/. The XgNT^X coupling is 
also desirable to generate acceptable fi and /x| terms (see section 5). 

In section it was also demonstrated that with a total of three mes- 
senger doublets the fine tuning required in electroweak symmetry breaking 
could be alleviated. By inserting the vev for $2 into W4, the doublet mass 
matrix is given as 



i^h 5g, Xg) 



/Ml 


\ Xsv^^ 



\ 

Xsv^s,^ 

/ 



( ^ \ 

\Xg J 



(65) 



and all other masses are zero. At Mqut the heavy doublets are (X/,X), 
(5g,Xg) and (5g,Xg), leaving the four vector-like doublets in (5/i,5/i), (5i,5i), 
(X,X;) and {Xh,Xh) massless at this scale. Of these four pairs, (5/i,5/j) are 
the usual Higgs doublets and the other three pairs couple to S. 

The (renormalizable) superpotential at scales below Mqut is then 



W = XgNq2q2 + -XnN^ + XioN'q3q2 

+XnN'q2q3 - XhNHuH^ + ^Xn'N'' 
+XqSIiIi + XiSqiqi + XSSI2I2 
+X,Sl3l3 + ^XsS' + Wr, 



(66) 



where the fields have been relabeled to make, in an obvious notation, their 
SU{3) X SU{2) X f/(l) quantum numbers apparent. 

We conclude this section with comments about both the choice of Z4 as 
a discrete symmetry and about non-renormalizable operators in our model. 

The usual i?-parity violating operators IOsm^sm^sm are not allowed by 
the discrete symmetries, even at the non-renormalizable level. In fact, im- 
parity is a good symmetry of the effective theory below Mqut- By inspection, 
the fields that acquire vevs at Mqut are either invariant under Z4 or have 
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a Z4 charge of 2 (for example, $1), so that a Z2 symmetry is left unbroken. 
In fact, the vevs of the other fields 5, N , N' and the Higgs doublets do not 
break this Z2 either. By inspecting the Z4 charges of the SM fields, we see 
that the unbroken Z2 is none other than the usual i?-parity. So at Mqut, the 
discrete symmetry Z4 is broken to Rp. We also note that the Z4 symmetry is 
suffucient to maintain, to all orders in 1/Mpi operators, the vev structure of 
$1 and $2, i-G-, to forbid unwanted couplings between $1 and $2 that might 
destabilize the vev structure [0]. This pattern of vevs was essential to solve 
the doublet-triplet splitting problem. It is interesting that both i?-parity 
and requiring a viable solution to the doublet-triplet splitting problem can 
be accommodated by the same Z4 symmetry. 

The non-SM matter fields {i.e., the messenger 5's and X's and the light 
triplets ) have the opposite charge to the SM matter fields under the unbroken 
Z2. Thus, there is no mass mixing between the SM and the non-SM matter 
fields. 

Dangerous proton decay operators are forbidden in this model by the 
discrete symmetries. Some higher dimension operators that lead to proton 
decay are allowed, but are suffuciently suppressed. We discuss these below. 

Renormalizable operators such as lOsAilOsM^g and lOsM^sM^q are for- 
bidden by the Z3 symmetries. This is necessary to avoid a large proton decay 
rate. A dimension-6 proton decay operator is obtained by integrating out the 
colored triplet scalar components of 5g or 5^. Since the colored scalars in 5g 
and 5q have a mass ~O(50 TeV), the presence of these operators would have 
led to an unacceptably large proton decay rate. 

The operators 10sm10sm10sm5sm/Mpi and IOsmIOsmIOsaSsm 
{^^/Mpi)"'/Mpi, which give dimension-5 proton decay operators, are also 
forbidden by the two Z3 symmetries. The allowed non-renormalizable op- 
erators that generate dimension-5 proton decay operators are suffuciently 
suppressed. The operator IOsm^^sm^Osm^smN' / {MpiY , for example, is al- 
lowed by the discrete symmetries, but the proton decay rate is safe since 
fAT/ ~ 1 TeV. 

The operators 10j5j$i(X or Xq)/Mpi could, in principle, also lead to a 
large proton decay rate. Setting $1 to its vev, the superpotential couplings, 
for example, Aij(f/fDjX(3) + QiLjX{3)) are generated with Xij suppressed 
only by Vis,j^/Mpi. In this model the colored triplet (scalar) components of X 
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and Xq have a mass rriq ~ 500 GcV, giving a potentially large proton decay 
rate. But, in this model these operators are forbidden by the discrete sym- 
metries. The operator lOibj^iX S/Mpi is allowed giving a four SM fermion 
proton decay operator with coefficient ~ v s / Mpi)^ / ~ lO^'^'^GeV^^. 
This is smaller than the coefficient generated by exchange of the heavy gauge 
bosons of mass Mqut-, which is ~ 9gut/^gut ~ 1/2 10~^^GeV"'^ and so this 
operator leads to proton decay at a tolerable rate. 

With our set of discrete symmetries, some of the messenger states and the 
light color triplets are stable at the renormalizable level. Non-renormalizable 
operators lead to decay lifetime for some of these particles of more than 
about 100 seconds. This is a problem from the viewpoint of cosmology, since 
these particles decay after Big-Bang Nucleosynthesis (BBN). With a non- 
universal choice of discrete symmetries, it might be possible to make these 
particles decay before BBN through either small renormahzable couphngs to 
the third generation (so that the constraints from proton decay and FCNC 
are avoided) or non-renormalizable operators. This is, however, beyond the 
scope of this paper. 

8 Conclusions 

We have quantified the fine tuning required in models of low energy gauge- 
mediated SUSY breaking to obtain the correct Z mass. We showed that 
the minimal model requires a fine tuning of order ~ 7% if LEP2 does not 
discover a right-handed slepton. We discussed how models with more mes- 
senger doublets than triplets can improve the fine tuning. In particular, a 
model with a messenger field particle content of three {I + lys and only one 
(g -|- q) was tuned to ~ 40%. We found that it was necessary to introduce 
an extra singlet to give mass to some color triplets (close to the weak scale) 
which are required to maintain gauge coupling unification. We also discussed 
how the vev and F-component of this singlet could be used to generate the 
/J, and B/j, terms. We found that for some region of the parameter space this 
model requires ~ 25% tuning and have shown that limits from LEP do not 
constrain the parameter space. This is in contrast to an NMSSM with one 
(5 -|- 5) messenger fields, for which we found that a fine tuning of ~ 1% is 
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Table 2: 5'[/(5) x ^t/ (5)' x Z3 x Z3 x Z4 quantum numbers for the fields of 
the model discussed in section 7. (a, 6, c) are the generators of Z3 x Z3 x Z4. 
The three SM generations are labeled by the index i. 
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required and that limits from LEP do significantly constrain the parameter 
space. 

We further discussed how the model with split messenger field represen- 
tations could be the low energy theory of a SU (5) x SU (5) GUT. A mecha- 
nism similar to the one used to solve the usual Higgs doublet-triplet splitting 
problem was used to split the messenger field representations. All operators 
consistent with gauge and discrete symmetries were allowed. In this model 
i?-parity is the unbroken subgroup of one of the discrete symmetry groups. 
Non-renormalizable operators involving non-SM fields lead to proton decay, 
but at a safe level. 
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10 Appendix 

In this section the Barbieri-Giudice parameters for both the MSSM and 
NMSSM in a gauge mediated SUSY breaking scenario are presented. 

In an MSSM with gauge mediated SUSY breaking, the fundamental pa- 
rameters of the theory (in the visible sector) are: A^ess] \] l-i'] a-nd /i^. Once 
electroweak symmetry breaking occurs, the extremization conditions deter- 
mine both m| and tan /3 as a function of these parameters. To measure the 
sensitivity of m| to one of the fundamental parameters Aj, we compute the 
variation in m| induced by a small change in one of the Aj. The quantity 



5m| 



c(m|; Xi) 



5\ 
Xi 



(67) 
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where 

measures this sensitivity 0. In the case of gauge mediated SUSY breaking 
models, there are four functions c(m|; Aj) to be computed. They are: 



2 2n _ 2/iV tan2/3 + l Atan^ (3{fll - fil) 

c[mz; /i ) - + (^^^2 ^ _ (^2 _ ~^l){t^^2 ;5 + 1) _ ^|(tan2 /3 - 1)) , 



(69) 



'^"^^'^^"^^ = ^tan^/j-l)3 
4 /i? - ^2 



tan^ /? 



for large tan/3, (70) 



, , , , ^ A? dmi dm% 



4 2 tan2/5 ^m^^^ / o'"?~'"2 tan^/J + l 



m| *tan2/3-l ^A? ' "/i? + /^i (tan^ /5 - l)s 
4 9m?r 



for large tan/5. (71) 



This measures the sensitivity of m| to the electroweak scale value of At, 
\t{Mweak)- The Yukawa coupling \t{Mweak) is not, however, a fundamental 
parameter of the theory. The fundamental parameter is the value of the 
coupling at the cutoff A° = Mqut or Mpi of the theory. We really should be 
computing the sensitivity of m| to this value of Aj. The measure of sensitivity 
is then correctly given by 

(73) 

We remark that for the model discussed in the text with three l+l and one g+ 
g messenger fields, the numerical value of (At(A°)/At(M^eafc))c^At(M^eafc)/f^At(A°) 
is typically ~ 0.1 because \t{Mweak) is attracted to its infra-red fixed point. 
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This results in a smaller value for c{'m\\ At) than is obtained in the absence 
of these considerations. 

With the assumption that rn^Hu ^"^^ "^^d scale with A^^^^, we get 

c("^I;AL.J = c(m|; m^J + c(m|; m^J 

tan^ 13+ I 
= 1+2^ ' 



ml (tan2/3-l)2 
4tan2/3(m|^^ + m|^J(/22 - jil)/ml 



(74) 



(yUi - 1^1) (tan^ (3 + 1)- mKtan^ (3-1)' 
The Barbieri-Giudice functions for mt are similarly computed. They are 

c{mt] nl) = ^c{ml; fij) + ^_^^^2p ^ (^5) 

c(m, ^^') = \c{ml, ^') + ~^,,Jf^ _ (76) 

c(-. AL J - ^c(.|; AL J - (78) 

Since niz and are measured, two of the four fundamental parameters 
may be eliminated. This leaves two free parameters, which for conveinence 
are chosen to be A^ess and tan (3. 

In a NMSSM with gauge mediated SUSY breaking, the scalar potential for 
iV, Hu and at the weak scale is specified by the following six parameters: 
Aj = m%,mjj^,mjj^, the NHuHd coupling A//, the scalar NHuHd coupling 
Ah, and the A^^ coupling, Atv- In minimal gauge mediated SUSY breaking, 
the trilinear soft SUSY breaking term NHuH^ is zero at tree level and is 
generated at one loop by wino and bino exchange. In this case, AniXi) = 
XHA{Xi). Since the trilinear scalar term A^'^ is generated at two loops, it 
is small and is neglected. The extremization conditions which determine 



= 5'|'i^^/4 (f = y vl + v'^) , tan P = Vu/vd and vn as a function of these 
parameters are given in section ^ Equation ^ can be written, using n = 
XrVn/V^ as 

A^ 111 

ml + 2-ff/ - XhXn-v^ sin 2/3 + -A>2 _ Ahv'^Xh sin 2/3 = 0. (79) 
Xjj 2 2 Afj, 
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Equation E3| is 



-22,2 2 



tan^ (3 



1 



0. 



(80) 



1 - tan2 (3 ^n- tan2 (3 
Substituting v% from equation in equation ^ and then using this expres- 
sion for /i| in equation |^ gives 

/ 9 9 9n ^ -^i? / 9 1 , 9 9n X / 2u 1 A?/ SiU 2/5 , 

(m^„+m^^+2/i2)sin2/?+-^(m^ + -A>2)+AH(-T^--^4 ^) = 0. 

(81) 

The quantity c = {\i / m%){dm% / d\i) measures the sensitivity of mz to these 



parameters. This can be computed by differentiating equations 79, 80 and 81 



with respect to these parameters to obtain, after some algebra, the following 
set of linear equations: 



{A + Aa^)X' = B' + B 
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(89) 
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(93) 



In deriving these equations AH^Xi) — AijA(Aj) was assumed and dA/dXn was 
neglected. Inverting these set of equations gives the c functions. We note 
that these expressions for the various c functions are vahd for any NMSSM 
in which the scalar term is negligible and the NHuH^ scalar term is 
proportional to Xh- In general, these 6 parameters might, in turn, depend on 
some fundamental parameters, Aj. Then, the sensitivity to these fundamental 
parameters is: 



Aj (9m| 



A. 



nir 



^ dXj drri^z 



, Xj ^ ^' '^dX, 



(94) 
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For example, in the NMSSM of section ^, the fundamental parameters are 
Amess, ^H, ^N, K and Ag {Afj is a function of Xh and A^ess)- Fixing mz and 
rrit leaves 3 free parameters, which we choose to be Amess, Xh and tan/?. As 
explained in that section, the effect of Xh in the RG scaling of m|^^ and rn'jj^ 
was neglected, whereas the sensitivity of to Xh could be non- negligible. 
Thus, we have 

c[mz] Xh) = c[mz] Xh) + c^m^; m^) . (95) 

We find, in our model, that c(m|; m^) is smaller than c(m|; A//) by a factor 
of ~ 2. Also, using approximate analytic and also numerical solutions to the 
RG equation for m^, we find that (A///m|r)(5m^/9A//) is ^ 0.1. Conse- 
quently, in the analysis of section || the additional contribution to c(m|; A^) 
due to the dependence of on Xh was neglected. A similar conclusion is 
true for Atv- Also, 

g(m|; Ag) = c(m|; m\) \ . (96) 

We find that {Xq/m\){dm\ / dXq) is ~ 1 so that c{m\;Xq) is smaller than 
c{w?z] A//) by a factor of 2. 
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Captions: 

Figure 1: Contours of c(m|; fi^) =(10, 15, 20, 25, 40, 60) for a MSSM with a 
messenger particle content of one (5 + 5). In figures (a) and (c) sgn{fi) = — 1 
and in figures (6) and (d) sgn{fi) = +1. The constraints considered are: (I) 
=75 GeV , (II) m^o + m^o = 160 GeV, (III) m^^ =85 GeV, and (IV) 
m^o + m^o = 180 GeV. A central value of rritop =175 GeV is assumed. 

Figure 2: Contours of c(m|; /x^) =(1, 2, 3, 5, 7, 10) for a MSSM with a 
messenger particle content of three {l + Jys and one (q + q). In figures (a) and 
(c) sgn{fi) = —1 and in figures (6) and (d) sgn{ii) = +1. The constraints 
considered are: (I) mg^ =75 GeV , (II) m^o + m^o = 160 GeV, (III) mg^ =85 
GeV, and (IV) m^o + m^^o = 180 GeV. A central value of rritop =175 GeV is 
assumed. 

Figure 3: Contours of c(m|; \h) for the NMSSM of Section 5 and a 
messenger particle content of three {I + lys and one {q-\- q). In figures (a) 
and (fe), c(m|; Aij)=(4, 5, 6, 10, 15) and \h =0.1. In figures (c) and (d), 
c(m|; Aij) =(3, 4, 5, 10, 15, 20) and Aij=0.5. The constraints considered 
are: (I) rrih + rUa — mz, (II) =75 GeV, (III) m^o + m^o — 160 GeV, 
(IV) ruh = 92 GeV, (V) mg^ =85 GeV, and (VI) m^o+m^o = 180 GeV. For 
Xh =0.5, the limit 70 GeV constrains tan/3^ 5 (independent of A^ess) 
and is thus not shown. A central value of rritop =175 GeV is assumed. 

Figure 4: Contours of c(m|; A^^) =(50, 80, 100, 150, 200) for the NMSSM 
of Section 5 with Xh =0.1 and a messenger particle content of one (5 + 5). 
The constraints considered are: (I) rrih + rria = mz, (II) ^ng^ =75 GeV, (III) 
m^o + m^o = 160 GeV, (IV) rrih =92 GeV, (V) mg^ =85 GeV, and (VI) 
m^o + m^o = 180 GeV. A central value of rritop =175 GeV is assumed. 
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